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ABSTRACT 

Two-dimensional  free  surface  flows  without  waves,  produced  by  a  submerged 
sink  in  a  reservoir,  are  computed  numerically  for  various  configurations.  For 
a  sink  above  the  horizontal  bottom  of  a  layer  of  fluid,  there  are  solutions 
for  all  values  of  the  Froude  number  F  greater  than  some  particular  value. 
However,  when  the  fluid  is  sufficiently  deep,  there  is  an  additional  solution 
for  one  special  value  of  F  <  1.  The  results  for  a  sink  at  the  vertex  of  a 
sloping  bottom,  treated  by  Craya  and  by  Hocking,  and  for  a  sink  in  fluid  of 
infinite  depth,  treated  by  Tuck  and  Vanden-Broeck ,  are  confirmed  and  extended. 
In  particular  it  is  shown  that  as  the  bottom  tends  to  the  horizontal,  the 
solution  for  a  sink  at  the  vertex  of  a  sloping  bottom  approaches  a  solution 
for  a  horizontal  bottom  with  F  =  1.  However  solutions  are  found  for  all 
values  of  the  Froude  number  F  >  1  for  a  sink  on  a  horizontal  bottom. 
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SIGNIFICANCE  AND  EXPLANATION 


When  fluid  is  withdrawn  from  a  reservoir,  the  free  surface  may  be  drawn 
down  as  in  Figures  1  and  6.  In  order  to  investigate  this  phenomenon  we 
consider  two-dimensional  free  surface  flows  without  waves,  produced  by  a 
submerged  sink  in  a  reservoir.  Numerical  solutions  are  obtained  for  various 
configurations.  For  a  sink  above  the  horizontal  bottom  of  a  layer  of  fluid, 
there  are  solutions  for  all  values  of  the  Froude  number  F  greater  than  some 
particular  value.  However,  when  the  fluid  is  sufficiently  deep,  there  is  an 
additional  solution  for  one  special  value  of  F  <  1.  We  were  led  to  look  for 
these  solutions  by  our  experience  with  other  free  surface  flows  with  gravity, 
such  as  flows  over  weirs  in  channels  and  flows  around  lips  of  teapot  spouts. 

In  those  cases  we  found  that  in  fluids  of  infinite  depth  there  was  a  flow  only 
for  a  special  value  of  the  appropriate  Froude  number.  This  kind  of  flow  also 
occurred  in  fluids  of  finite  depth,  but  in  addition  there  were  solutions  for 
all  Froude  numbers  greater  than  some  particular  value.  €)ur  present  results 
show  that  this  is  also  the  case  for  free  surface  flows  produced  by  sinks. 


The  reponsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


FREE  SURFACE  FLOW  DUE  TO  A  SINK 


*  1  ^ 
Jean-Marc  Vanden-Broeck  '  and  Joseph  B.  Keller 


‘,2 


§1.  Introduction. 

When  fluid  is  withdrawn  from  a  reservoir  by  a  sink  of  strength  Q  at  depth  ft,  the  surface 
above  the  sink  may  be  drawn  down,  as  in  Figures  1  and  6.  An  exact  solution  of  this  type 
was  found  by  Craya  (1949)  (also  in  Yih  (1965),  pp.  124-126)  when  the  bottom  sloped 
downward  at  the  angle  =  7r/3  from  the  vertical,  and  a  numerical  solution  was  found 
by  Tuck  and  Vanden-Broeck  (1984)  for  ^  =  0.  Then  Hocking  (1985)  obtciined  numerical 
solutions  for  a  sequence  of  angles  ranging  from  0  to  7r/2.  In  each  case  there  was  just  one 
solution  without  waves.  We  have  recomputed  these  two-dimensional  flows  and  confirmed 
their  solutions. 

For  a  horizontal  bottom,  however,  we  have  found  solutions  for  all  values  of  F  greater 
than  some  particular  value.  When  the  fluid  is  sufficiently  deep,  there  is  an  additonal 
solution  for  one  special  value  of  F"  <  1.  We  were  led  to  look  for  these  solutions  by  our 
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experience  with  other  free  surface  flows  with  gravity,  such  as  flows  over  weirs  in  channels 
and  flows  around  lips  of  teapot  spouts.  In  those  cases  we  found  that  in  fluid  of  infinite 
depth  there  was  a  flow  only  for  a  special  value  of  the  appropriate  Froude  number.  This 
kind  of  flow  also  occurred  in  fluid  of  finite  depth,  but  in  addition  there  were  solutions  for 
all  Froude  numbers  greater  than  some  particular  value.  Our  present  results  show  that  this 
is  also  the  case  for  free  surface  flows  produced  by  sinks. 

In  section  2  we  compute  wave-free  flows  produced  by  a  sink  in  a  liquid  layer  of  finite 
depth  with  a  horizontal  bottom.  We  find  that  the  sink  strength  Q  and  the  Froude  number 
F  are  determined  by  the  sink  depth  h  and  the  distance  W  from  the  sink  to  the  bottom. 
The  result  are  discussed  in  section  3.  It  is  pointed  out  that  the  solutions  are  subcritical 
for  0  <  h/W  <  0.70  and  supercritical  for  0.70  <  h/W  <  0.76,  while  no  solutions  are  found 
for  h/W  >  0.76.  Additional  supercritical  solutions  are  constructed  in  section  4,  and  it  is 
indicated  that  there  are  additional  subcritical  solutions  with  waves.  In  section  5  we  find 
solutions  for  a  sink  at  the  corner  on  a  sloping  bottom,  which  were  mentioned  in  the  first 
paragraph  of  this  Introduction. 


§2.  Sink  above  a  horizontal  bottom. 

Let  us  consider  the  flow  in  the  region  of  the  2-plane  shown  in  Figure  1.  At  the  distance  W 
from  a  horizontal  bottom  there  is  a  sink  with  strength  2Q.  The  bottom  BI  is  a  streamline 
on  which  we  require  that  the  stream  function  tp(xy  y)  =  Q.  The  portion  DS  of  the  vertical 
wall,  where  i  =  0,  0  <  y  <  W,  is  part  of  the  same  streamline.  The  portion  SC  of  the 
verticad  wall,  where  x  =  0,  W  <  y  <  yc,  is  part  of  another  streamline  on  which  t/>(x,  y)  =  0. 
That  streamline  continues  along  the  free  surface  Cl,  which  leaves  the  wall  tangentially  at 
some  point  C  which  is  to  be  determined. 

We  denote  by  h  the  depth  of  the  sink  below  the  level  of  the  free  surface  at  infinity. 
We  also  choose  the  unit  of  length  and  the  unit  of  velocity  so  that  Q  =  1  and  y  =  1.  Then 
we  introduce  the  complex  velocity  potential  f{z)  =  (p{x,y)  iil'(x,y).  In  the  /-plane,  the 
flow  r^'gion  is  the  strip  0  <  V’  <  1  with  the  streamline  ICS  on  j/’  =  0  and  the  streamline 
IBS  on  0  =  1.  We  map  it  onto  the  lower  half  of  the  unit  circle  of  the  auxiliary  t-plane  by 
the  transformation 


/=  -llog 

I 


4t 


(2.1) 


(See  Figure  2). 
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Next  we  consider  the  complex  velocity  ^(t)  =  u  —  tv  as  a  function  of  t  in  the  semicircle 
in  Figure  2.  Here  u  =  and  v  =  ipy  are  the  velocity  components.  Since  there  is  a  sink  of 
strength  -2  at  <  =  0  and  a  right  angled  corner  at  x  =  y  =  0,  C  niust  have  singularities  at 
these  two  points  and  be  regular  elsewhere.  The  appropriate  singularities  are  proportional 
to  at  <  =  0  and  to  {t  4-  b)2  at  <  =  -b,  where  -b  is  the  image  of  the  corner  B  in  the 
t-plane  (see  Figure  2).  Therefore  we  write  C(0  in  the  form 

1  °° 

C(t)  =  -i(t  +  b)2f“^  exp(^a„f").  (2.2) 

n=0 


The  coefficients  a„  are  to  be  found. 

On  the  wall  from  B  to  C,  u  =  0  so  )  must  be  imaginary  for  —  6  <  <  <  1.  In  addition 
on  the  bottom  from  B  to  /,  u  =  0  so  <^{t)  must  be  real  for  —1  <  t  <  —b.  Both  of  these 
conditions  are  satisfied  by  (2.2)  if  we  choose  all  the  a„  to  be  real. 

On  the  free  surface  from  C  to  I  the  pressure  is  assumed  to  be  constant.  By  using  the 
Bernoulli  equation  we  can  write  this  condition  as  |^l^+2y  =  constant.  In  the  <-plane  it  must 
hold  on  the  circular  arc  t  =  e“^,  0  >  <7  >  -n.  Differentiating  this  condition  with  respect  to 
<7  yields  d\C,\'^/d(7  +  25y/5(7  =  0.  Now  on  the  circular  arc  (2.1)  yields  /  =  — ^  logcos((7/2), 
which  is  real,  so  f  =  y?  there.  Differentiating  this  relation  gives  d^pjdc  =  ^tan(<7/2). 
We  can  use  this  expression  to  write  dy/da  —  (dy  f  d<p){dip / dcr)  =  [v/(u^  +  w  tan(a/2). 
Upon  inserting  this  result  into  the  differentiated  form  of  the  constant  pressure  condition, 
we  obtain 


W  , 

da 


V 


=  0 


on 


t  =  e 


i<T 


0  >  a  >  — TT. 


(2.3) 


We  now  set  t  =  in  (2.2)  to  get  and  we  substitute  that  expression  into  (2.3). 

WV  will  use  the  resulting  equation  to  determine  the  unknown  coefficients  a„  that  occur  in 

(2.2). 

To  do  so,  we  truncate  the  infinite  .series  in  (2.2)  after  JV  —  1  terms.  We  find  the  JV  —  1 
coefficients  a„  and  the  constant  b  by  collocation.  Thus  we  introduce  iV  —  1  mesh  points. 
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<7/ = -7r(7-  i)/(.V-l),  1  =  1 . .V-1.  (2.4) 

By  using  (2.2)  we  obtain  C  and  0\C,\^ /Oa  at  a  =  a/  in  terms  of  the  coefficients  a„  and 
f).  Uinm  substituting  these  expressions  into  (2.3)  we  obtain  N  -  1  nonlinear  algebraic 


4 


equations  for  the  N  unknowns  and  b.  An  iV-th  equation  is  obtained  by  fixing 

the  position  of  the  sink.  This  is  achieved  by  specifying  either  a  value  for  b  or  a  value  for 
h/W.  To  fix  h/W  we  evaluate  h  and  W  by  integrating  the  identity 

dx  dy  1 

along  the  real  diameter  and  along  the  circle  in  the  t-plane. 

We  solved  the  resulting  system  of  N  equations  by  Newton’s  method.  Once  this  system 
was  solved  we  obtained  the  shape  of  the  free  surface  by  integrating  (2.5)  numerically.  In 
addition,  we  calculated  the  Froude  number  F  based  on  the  total  depth  at  infinity,  which 
is  defined  by 

F  =  Q[g(H  +  W)T'^.  (2.6) 


From  our  choice  of  dimensionless  variables,  it  follows  that  F  is  given  by 

F  =  [C(-l)]l  (2.7) 

The  coefficients  were  found  to  decrease  rapidly.  For  example,  a\  ^  0.3,  aio  ~ 
—2  X  10“®  and  020  ~  ~2  x  10~*  for  h/W  =  0.74.  Most  of  the  calculations  were  done  with 
N  =  30. 


§3.  Discussion  of  results. 

The  present  problem  is  qualitatively  similar  to  the  flow  over  a  weir  which  we  have  alu  ady 
studied  (Vanden-Broeck  &  Keller  (1986)).  In  particular,  the  flux  Q  and  the  ratio  hjW 
cannot  be  specified  independently  when  h/W  is  sufficiently  small. 

Following  our  analysis  of  weir  flows  we  define  a  ^'discharge  coefficient"  Cih/W]  by 
the  relation 

Q  =  C{h/W)g'^hh  (3.1) 

\^e  assume  that  for  h/W  small,  C(h/W)  can  be  expanded  in  a  finite  Taylor  scries: 

C(h/W)  =  C(0)  +  ^C'(Q)  +  0[(h/W)^]  (3.2) 

Upon  substituting  (3.1)  into  (2.6)  and  expanding  for  h/W  small,  we  obtain 

3 

iciO)  +  -  IC(0)]  +  0[{h/W)^]]  . 


(3.3) 


Relation  (3.3)  is  similar  to  the  relation  (2.6)  we  derived  for  weir  flows  (Vanden-Broeck  & 
Keller  (1986)). 

Tuck  ^md  Vanden-Broeck  (1984)  solved  the  sink  problem  numerically  for  h/W  = 
0.  They  found  a  unique  solution  characterized  by  {2Q)'^{gh^)~^  =  12.622.  This  result, 
together  with  (3.1),  implies 

C(0)  =  1.776.  (3.4) 

In  order  to  find  C'(0)  and  C{-^)  for  h/W  0  we  plotted  our  numerical  results  for 
F{h/W)~^  versus  h/W.  For  h/W  small  the  graph  is  close  to  a  straight  line  of  slope  —2.4 
which  intersects  the  F(h/W)~^^^  axis  at  1.78.  Therefore  (3.3)  implies  that 

C(0)  ~  1.78  (3.5) 


Relations  (3.5)  and  (3.6)  yield 


C'(0)  -  |C(0)  ~  -2.4 


C'(0)  =  0.3. 


We  note  that  (3.5)  agrees  with  the  value  (3.4)  previously  obtained  by  Tuck  and  Vanden- 
Broeck  (1984). 

In  Figure  3  we  have  plotted  the  values  of  F  versus  h/W  and  versus  b.  For  h/W  <  0.70, 
which  corresponds  to  6  >  0.37,  the  flow  is  subcritical,  t.e.,  F  <  1.  For  h/W  >  0.70, 
corresponding  to  6  <  0.37,  the  flow  is  supercritical,  ».e.,  F  >  1.  As  h/W  — f  0.76,  F  00. 
These  solutions  exist  only  for  h/W  <  0.76. 

As  h/W  increases  from  0  to  0.76,  the  parameter  b  decreases  monotonically  from  1  to 
0.33.  We  shall  show  later  that  this  value  is  exactly  1/3.  A  typical  free  surface  profile  for 
h/W  =  0.5  is  shown  in  Figure  1. 

Our  subcritical  solutions  for  h/W  <  0.70  are  characterized  by  a  uniform  stream  at 
infinity.  This  is  exceptioned  since  usually  a  subcritical  flow  has  a  train  of  waves  at  infinity. 
We  believe  that  subcritical  solutions  with  waves  also  exist  for  the  present  problem.  Their 
wave  amplitude  a  must  be  a  function  of  F  and  h/W : 

a  =  a{F,h/W).  (3.8) 

The  special  expression  (2.2)  requires  that  there  be  no  waves  at  infinity.  Relation  (3.8) 
indicates  that  this  is  possible  only  if  some  relation  between  F  and  h/W  is  satisfied,  namely 
a(F,  h/W)  =  0.  Such  a  relation  is  exactly  what  we  discovered  numerically  (see  Figure  3). 
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Tlie  precedin.;:  considerations  suggest  that  there  are  also  additional  supercritical  so- 
'.uiiuns  for  h/W  >  0.70 .  Supercritical  flows  are  characterized  by  the  presence  of  exponen¬ 
tially  decreasing  terms  at  infinity.  Thus  generally  ^  has  the  form 

~  (oo  +  as  -oo.  (3.9) 


Here  X  is  the  smallest  positive  root  of 


ttA  —  F  ^tan7rA  =  0. 


(3.10) 


■ ^  ■ 


Equation  (3.10)  has  real  solutions  only  for  F  >  1,  t.e.,  only  for  supercritical  flow.  For 
F  <  1,  (3.10)  has  purely  imaginary  roots  and  then  (3.9)  corresponds  to  a  train  of  small 
amplitude  waves  at  infinity.  Therefore  the  exponentially  small  term  which  occurs  in  the 
expression  for  supercritical  flows  is  the  “analytic  continuation''  of  the  small  amplitude  wave 
term  of  subcritical  flows.  As  in  the  subcritical  case,  the  factor  A  in  (3.10)  must  depend 
upon  F  and  h/W: 

A  =  A(F,h/W).  (3.11) 

In  terms  of  the  transformation  (2.1)  from  /  to  t,  we  can  rewrite  (3.9)  as 

C  ~  Coo  +  A(<  4- 1)"^,  t  — 1.  (3.12) 

The  special  expression  (2.2)  does  not  allow  a  singularity  like  that  in  (3.12)  as  t  —1. 
Therefore,  the  supercritical  solutions  we  have  computed  correspond  to  A(F,  h/W)  =  0. 
This  agrees  with  our  numerical  results  which  show  that  solutions  of  the  form  (2.2)  can  be 
obtained  only  if  some  relation  between  F  and  h/W  is  satisfied,  as  Figure  3  shows. 

§4.  Additional  supercritical  flows. 

We  shall  now  construct  solutions  for  which  A{F,h/W)  ^  0.  For  F  =  oo,  the  velocity  is 
constant  on  the  free  surface  and  the  problem  has  an  exact  solution,  namely 

=  (4,1) 

t(  1  t6)  2 

Solution  (4.1)  is  defined  for  all  values  of  6.  By  using  (2.1)  and  (4.1)  we  find  that 

irtp  _  36 

c  ~  C(-l)  +  e  2  ^-((-1)  as  0-^-oo.  (4.2) 


V-V  M 


r  ^  > 
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s.*. 


Relation  (3.10)  shows  that  A  — »  |  as  F  — »  oo.  Therefore  matching  (3.9)  and  (4.2),  and 
using  (2.7),  yields 

2  1  —  OL 

A F3  - — —  as  F oo.  (4.3) 

2  —  2b 

The  asymptotic  expression  (4.3)  vanishes  for  b  =  1/3.  The  corresponding  value  of  Ji/W  is 
0.76.  This  is  the  solution  for  F  =  oo  obtained  by  using  (2.2).  (See  Figure  3).  However 
solutions  for  F  =  oo  exist  for  b  ^  1/3.  They  are  given  by  (4.1). 

Solutions  for  F  large  can  be  obtained  by  using  perturbation  theory:  relation  (4.1)  can 
be  considered  as  the  first  term  of  an  expansion  in  powers  of  F~^ .  Instead  of  doing  that  we 
shall  solve  the  problem  numerically  for  arbitrary  values  of  F. 

We  generalize  (2.2)  to  allow  the  singularity  (3.12)  at  t  =  —1.  Therefore  we  write 


c  = 


1  , 

—i(t  +  b)'it  ^  exp 


j4(1  +  t)^'^  ^  ^  Unt" 

n=0 


(4.4) 


Here  A  is  the  smallest  positive  root  of  (3.10).  The  coefficients  a„  and  the  constant  A  are 
to  be  found.  We  truncate  the  infinite  series  after  N  —  1  terms  and,  as  before,  satisfy  (2,3) 
at  the  collocation  points  (2.4).  This  leads  to  N  —  1  algebraic  equations  for  the  A'^  +  1 
unknowns  6  and  A.  Two  more  equations  are  obtained  by  specifying  F  and  the 

position  of  the  sink,  i.e.,  h/W  or  b.  The  resulting  system  of  A^  +  1  equations  with  N  f-  1 
unknowns  is  solved  by  Newton’s  method. 

We  note  that  we  specify  two  parameters,  F  and  h/W  or  b,  in  order  to  obtain  u  unique 
solution.  We  checked  numerically  that  specifying  more  or  less  parameters  did  not  yield 
convergence  for  N  large. 

2 

In  Figure  4  we  present  numerical  values  of  the  parameter  A/F3  appearing  in  (4.1) 
versus  F  for  b  =  0.3,  0.35  and  0.4.  As  F  —*  oo,  the  curves  approach  asymptotically  tl;.: 
values  (4.3).  The  points  at  which  the  curves  intersect  the  F-axis  correspond  to  supercritical 
solutions  for  which  =  0  in  (4.4).  These  are  the  solutions  described  by  the  iiortion 
j  <  6  <  0.37  of  the  curve  in  Figure  3. 

For  6  <  1/3  or  5  >  0.37  the  curves  in  Figure  4  do  not  intersect  the  F-a.vis.  For 
h  >  0.37  we  found  that  as  the  Froude  number  is  decreased  from  infinity  each  ciir’-c  in 
Figure  4  reaches  a  limit  point  and  then  turns  back.  Therefore  for  some  valuc.s  0"  the 
Froude  number  two  different  solutions  are  possible.  A  typical  profile  for  h  —  0.4,  F  -  1.3 
and  AF~^f^  =  —1.15  is  shown  in  Figure  5.  We  expect  that  each  solution  branch  f<u 
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a  ;>  0.37  will  end  when  a  stagnation  point  appears  on  the  free  surface  with  a  120®  angle  at 
it. 

For  b  <  0.37  the  curves  of  Figure  4  extend  from  F  =  oo  to  F  =  1+.  We  did  not 
CO  -iii  ae  these  branches  into  the  range  F  <  1.  Since  A  ^  0  as  F  — ►  1+,  we  expect  (by 
anal^uc  continuationj  the  solutions  for  F  <  1  to  have  waves  at  infinity,  i.e.,  to  have  a  ^  0 
lu  (3.8). 

The  results  above  show  that  there  are  two  different  classes  of  supercritical  flows.  The 
first  class  is  chaxacterized  by  the  fact  that  subcritical  flows  Eu-e  ultimately  reached  when 
the  Froude  number  is  decreased  from  infinity  (curves  a  amd  b  in  Figure  4).  In  the  second 
class  a  limit  point  is  reached  before  the  subcritical  regime  is  reached  (curve  c  in  Figure  4). 
Another  example  of  supercritical  flow  of  the  second  class  htis  been  found  by  Vanden-Broeck 
(1986)  in  his  analysis  of  flows  under  a  gate. 

§5.  Sink  on  a  sloping  bottom. 

We  shall  now  extend  the  procedure  of  Section  2  to  obtain  the  flow  in  the  region  of  the 
z-plane  shown  in  Figure  6.  The  rigid  bottom  slopes  at  an  angle  0  from  the  vertical.  At 
the  origin  there  is  a  sink  of  strength  —2vQf{ir  —  13).  We  choose  rp  =  0  on  the  verticad  wail 
and  on  the  free  surface.  It  follows  from  the  value  of  the  sink  strength  that  ip  =  Q  on  the 
sloping  bottom.  As  before  we  choose  the  unit  of  length  and  the  unit  of  velocity  so  that 
Q  =  1  and  g  —  By  using  the  transformation  (2.1)  we  map  the  flow  domain  into  the 
interior  of  the  unit-circle  in  the  <-plane  (see  Figure  2). 

Since  there  is  a  sink  of  strength  —2nQf(ir  —  0)  at  t  =  0  and  a  source  at  t  =  —1, 
^  must  have  singularities  at  these  two  points  and  be  regular  elsewhere.  The  appropriate 
singularities  are  proportional  to  at  <  =  0  £ind  to  at  <  =  —1,  where  7  =  2  • 

Therefore  we  write 

00 

C  =  — 1(1  +  exp  •  (5.1) 

.n=0 

The  coefficients  a„  are  found  by  following  the  numericEd  procedure  of  Section  2. 

We  shall  consider  first  the  special  CEise  =  n/6  corresponding  to  7  =  1/3.  Then  by 
choosing  oo  =  3  log(^)  and  a„  =  0  for  n  >  1  we  obtain  an  exact  solution  of  (2.3).  This  is 
the  solution  found  by  Craya  (1949).  For  /3  =  0  (i.e.,  7  =  j)  the  problem  reduces  to  the 
flow  calculated  by  Tuck  and  Vanden-Broeck  (1984).  Hocking  (1985)  obtained  numerical 
solutions  for  a  sequence  of  angles  rsuiging  from  0  to  ^ .  Our  numerical  results  were  found 


Figure  6 


Vertical  section  of  a  reservoir  with  a  sloping  bottom  and 
a  vertical  wall.  The  sink  s  is  at  the  corner.  The 
free  surface  is  Cl ,  and  B  is  the  angle  between  the 
bottom  and  the  vertical.  The  figure  is  an  actual  profile 
for  B/tt  =  0.4.  The  vertical  scale  is  the  same  as  the 
horizontal  scale. 
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to  be  in  good  agreement  with  theirs.  In  particular  we  confirmed  that  a  unique  solution 
exists  for  each  value  oi  0  <  jS  <  All  these  solutions  axe  characterized  by  a  stagnation 
point  at  infinity.  For  0  =  ^  there  is  a  uniform  stream  at  infinity.  Our  numerical  results 
show  that  a  unique  solution  characterized  by  F'  =  1  is  obtained  as  yS  — +  ^ .  However  there 
is  a  family  of  solutions  for  (3  =  This  family  can  be  computed  by  choosing  b  =  0  in  (4.4) 
and  using  the  numerical  procedure  of  Section  2.  As  F  — ♦  1,  A  — >  0.  Therefore  the  limit 
as  F  — ♦  1  of  (4.4)  with  6  =  0,  agrees  with  the  limit  of  (5.1)  as  7  — »  1/2.  The  coefficients 
a„,  n  >  1  in  both  expressions  become  identical  and  the  difference  between  the  coefficients 
oo  is  equzd  to  A.  These  results  are  illustrated  in  Figures  7  and  8.  In  Figure  7  we  present 
numerical  values  of  the  depth  he  of  the  cusp  under  the  level  of  the  free  surface  at  infinity 
(see  Figure  6)  versus  F^  for  0  =  In  Figure  8  we  present  numeric2d  values  of  he  versus 
for  0  <  /3  <  f. 
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Figure  8.  Cusp  depth  h  versus  B/w  for  the  solutions  found  by 

Hocking.  Craya's  solution  is  at  6  =  —  and  that  of 
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Tuck  and  Vanden-Broeck  is  at  6  =  -z-  . 
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